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AIRCRAFT SYSTEM MODELING ERROR 
AND CONTROL ERROR 

ORIGIN OF THE INVENTION 

5 

This invention was made, in part, by one or more employ- 
ees of the U.S. government. The U.S. government has the 
right to make, use and/or sell the invention described herein 
without payment of compensation, including but not limited 
to payment of royalties. 10 

FIELD OF THE INVENTION 

This invention relates to minimization of modeling error 
and control error for an aircraft system. 15 

BACKGROUND OF THE INVENTION 

If an aircraft/spacecraft vehicle encounters a failure (such 
as a jammed control surface or loss of a surface), most con- 20 
trailers cannot adapt to the failure and a crash may occur. In 
most cases, the vehicle has enough redundant actuation 
mechanisms to salvage the vehicle. Several airplane crashes 
have occurred in the past where the pilot is unable to control 
the damaged airplane due to the pilot’s inability to learn to fly 25 
this altered aircraft configuration in the very short time avail- 
able. The flight computer, however, may have the necessary 
information as well as bandwidth available to learn the new 
dynamics, and control the vehicle within a reasonable time 
interval. 30 

The flight computer needs an intelligent controller that flies 
the vehicle with the baseline controller during nominal con- 
ditions, and adapts the design, when the vehicle suffers dam- 
age. Thus, given the information about the vehicle from all the 
available sensors, the control system needs to determine 35 
whether the vehicle is in its nominal state or is damaged. One 
approach to deal with this is to utilize smart algorithms that 
attempt to identify the vehicle characteristics and to change 
the control system, if necessary. This approach is known as 
Indirect Adaptive Control. For systems such as airplanes, 40 
there is usually very little time available to make changes to 
the control system, and this indirect approach is often insuf- 
ficient to achieve the desired safety metrics. Another 
approach, known as the direct adaptive control (“DAC”), 
looks directly at the errors, and updates the control law 45 
accordingly. This is typically much faster and meets the tim- 
ing considerations for airplane system implementations. 

The current state of the art implementation consists of the 
Intelligent Flight Control Architecture that uses a DAC 
approach. This has been implemented by us at the NASA 50 
Ames Research Center, and has been test flown on the F-15 
research aircraft at the Dryden Flight Research Center. The 
update law uses tracking error to change the control law. This 
approach is based on the work at the Georgia Tech Aerospace 
Engineering Department, under R. T. Rysdyk and A. J. Calise, 55 
“Fault Tolerant Flight Control Via Adaptive . . . Augmenta- 
tion” AIAA 98-4483. 

When operating in the real world, an airplane will always 
have tracking errors associated with its states. For example, 
when an pilot provides an aggressive stick command, there is 60 
always a large transient tracking error that eventually disap- 
pears as the controller continues to perform. Adaptation 
should typically occur only when the aircraft experiences 
damage or change in its flight configuration, which the origi- 
nal control design cannot deal with. Usually much effort goes 65 
into the design of the nominal baseline control design, which 
should be changed only if necessary. 


2 

What is needed is an approach that implements DAC that 
looks not just at the tracking error, but rather its characteris- 
tics or evolution over time to determine whether the controller 
needs to be adapted or left alone. The time evolution of the 
tracking error provides clues for investigating whether the 
system is in good health or has undergone damage/faults. This 
crucial piece of available information remains un-utilized in 
all the existing DAC approaches. 

SUMMARY OF THE INVENTION 

This invention presents a novel stable discrete-time adap- 
tive law that is designed and implemented for flight control to 
targets damages/modeling errors in a direct adaptive control 
(DAC) framework. The approach is based on the observation 
that, where modeling errors are not present, the original con- 
trol design has been tuned to achieve the desired performance. 
The adaptive control should, therefore, work towards achiev- 
ing back this performance only when the design has modeling 
uncertainties/errors or when the vehicle suffers damage or 
substantial flight configuration change. In this work, the base- 
line controller uses dynamic inversion with proportional-in- 
tegral augmentation. Dynamic inversion is carried out using 
the assumed system model. On-line adaptation of this control 
law is achieved by providing a parameterized augmentation 
signal to a dynamic inversion block. The parameters of this 
augmentation signal are updated to achieve the nominal 
desired error dynamics. Contrary to the typical Lyapunov- 
based adaptive approaches that guarantee only stability, the 
current approach investigates conditions for stability and for 
performance. A high-fidelity F-15 aircraft model is used to 
illustrate the overall approach. 

Operationally, the aircraft plant dynamics is modeled, 
using the original plant description without changes, and the 
parameters representing the plant components are monitored. 
Under normal conditions, the controller responds to an excur- 
sion in the tracking error e(k), which is the difference between 
the desired and the actual aircraft behavior, and drives this 
tracking error toward a zero value according to n asymptote 
curve that is characteristic of the controller. If the tracking 
error does not conform to, or lie close to, this asymptotic 
curve, a resulting error (difference between desired error 
behavior and actual error behavior) is observed. This differ- 
ence, called the performance error E(k), represents a differ- 
ence between normal aircraft parameters and damaged air- 
craft parameters, and its components are monitored. 

Assume that the system senses that (1) at least one com- 
ponent of aircraft tracking error e(k) is experiencing an excur- 
sion and (2) the return of this component value toward a 
reference value (e.g., a constant, such as 0) is not proceeding 
according to the expected controller characteristics (which 
gives rise to a non-zero magnitude IE(k)l above an expected 
threshold magnitude). Only when both conditions (1) and (2) 
are satisfied will the system reactivate the neural network 
(NN), change the plant dynamics according to the NN, and 
change the modeling of aircraft plant operation. Where con- 
dition (1) is satisfied but the return of the vector component 
e(k) toward the reference value proceeds according to the 
controller characteristics (E(k)=0), or within a selected 
neighborhood of this asymptote, so that condition (2) is unsat- 
isfied, the system will not change modeling of the plant opera- 
tion. In this latter instance, the NN will continue to model 
operation of the aircraft plant according to the original model . 
In a prior art approach, as long as condition (1) is satisfied, 
modeling of the aircraft plant dynamics is changed, irrespec- 
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tive of whether the components of the vector E(k) are follow- 
ing the controller characteristics. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

FIG. 1 illustrates the overall adaptive control architecture. 

FIGS. 2 A and 2B graphically illustrate pilot longitudinal 
and lateral stick input signals. 

FIGS. 3A-3C through 5A-5C graphically illustrate aircraft 
pitch, roll and yaw behavior for a selected test. to 

FIGS. 6A-6C and 7A-7D graphically illustrate perfor- 
mance errors for the roll, pitch and yaw behavior for a 
selected test. 

FIGS. 8A-8C through 10A-10C illustrate pitch, roll and 
yaw acceleration behavior, where a system failure is inserted 15 
at an intermediate time (10 sec). 

FIGS. 11A-11C illustrate roll, pitch and yaw performance 
errors corresponding to the system failure introduced in 
FIGS. 8A-C, 9A-C and 10A-C. 

FIGS. 12A-12D illustrate aircraft surface commands (sta- 20 
bilizer, rudder, canard) corresponding to the system failure 
introduced in FIGS. 8A-C, 9A-C and 10A-C. 

DESCRIPTION OF THE INVENTION 

25 

An adaptive controller, according to the invention, updates 
the nominal baseline control approach only if there is a mod- 
eling error or damage occurs or a substantial change in flight 
configuration occurs that cannot be corrected in a conven- 
tional manner by the controller. 30 

Control Architecture. 

FIG. 1 presents the overall adaptive control architecture 
and schematically illustrates an embodiment of adaptive con- 
trol architecture for practicing the invention. A reference 
module 1 1 provides a vector y r£?/ (k+ 1 ) of one or more aircraft 35 
system variables that are to be monitored and controlled, 
where k is a monotonically increasing time index. The refer- 
ence vector y re/ (k) is received at a first difference module 13 
that forms a difference vector, y re y(k)-y(k)=e(k) that repre- 
sents an error vector between a presently sensed aircraft sys- 40 
tern vector y(k) and the reference vector y re j(k), which 
becomes an input signal for a proportional integral augmen- 
tation (PI A) module 15. An output signal y Pe (k+l) the PI A 
module 15 is received at a first positive terminal of a second 
difference module 17, which receives the reference vector 45 
y re /k) at a second positive input terminal and receives an 
adaptive augmentation signal y^(k+l) at a negative input 
terminal. The output (difference) signal y^(k+l) of the sec- 
ond difference module 17 is received as an input signal at a 
dynamic inverse module 19, which generates an output con- 50 
trol signal vector u(k). The control signal vector u(k) is 
received by an aircraft physical plant module 21 that gener- 
ates an aircraft system output vector y(k+ 1 ), which is received 
by the first difference module 13 and as a first input signal y(k) 
at an adaptive augmentation (AA) module 23. The AA mod- 55 
ule 23 also receives a second input signal y Pe (k+l) and a third 
input signal y reJ (k+\) and generates an output adaptive aug- 
mentation signal y arf (k+l). 

The control system is given a command, y com (k+\) (e.g., 
pitch rate command from the pilot’s stick). The time index 60 
(k+ 1 ) refers to the desired value at the next time index (k+ 1 ). 
Given the knowledge of how fast or slow the aircraft plant can 
handle such a command, it is typically taken through a second 
order reference model, with appropriate damping and natural 
frequency to obtain the corresponding achievable reference 65 
command y re y(k+l). It is important to note that the value for 
the time index (k+1) for this reference signal is not necessar- 


4 

ily computed at time index (k+1), but is the desired reference 
value of the output the time index (k+1) that is computed at 
time index (k). The controller is designed to achieve a pre- 
scribed second order error dynamics with respect to this ref- 
erence command. Let this error dynamics, in a discrete form, 
be given in scalar form as: 


e(k+ l)+K Pe e(k)+K 1 e e { (k)=0 , 

(1A) 

e(k)=y(k)-y re /k) 

(IB) 


where e lk represents the integrated error until time index k. 
K Pe and K le are gains, chosen appropriately to have the 
desired transient response characteristics. Equation (1), with 
the definition of the error e(k), is used to compute the control 
input to achieve the desired error dynamics as follows. Equa- 
tion (1) can be re-expressed as 

y(k+l)=y re /k+l )+K Pe {y re /k)-y(k)}-K le e l (k) (2) 

The plant output y(k+ 1) must satisfy Eq. (2) to achieve the 
prescribed second order error dynamics. The right hand side 
of Eq. (2) can thus be labeled as y^(k+l), the desired plant 
output. Thus, 

^(*+ 1 )^ ( 3 ) 

Again, note that this value of the desired output at time index 
(k+1) is computed at time index (k). Let the plant dynamics be 
given as: 

y(k + 1) = f{y(k), y(k- 1) yik-py, u{k - 1), . . . ,u(k- p„)) ++ (4) 

g{y(k), y(k- 1), ... , y{k -p y \ u(k - 1), ... , u(k - p u )} u(k) 

We can thus invert the dynamics represented by Eq. (4) to 
compute the control function u(k) to achieve the desired error 
dynamics, Eq. (1), as: 


«(*) = {yrtfik + 1) + K P '{y„f(k) - y(k)) - K le edk) - < 5 ) 

y(k- 1), ... , y{k-p y \ u(k- 1), ... ,u(k-p u )}}/ 
g{y(k), y(k-l), ... ,y(k-p y ;u(k- 1), ... , u(k-p u )}, 

where f and g are functions characterizing the plant. 

This control input, with exact knowledge of the plant (f and 
g), will help achieve the desired second order error dynamics. 
With modeling uncertainties and other errors, we will not 
know f and g exactly, but only their estimates given by the 
model, F and g~. The adaptive augmentation is now designed 
to offset these modeling errors, so that we can get the same 
error dynamics or the desired performance. With the adaptive 
augmentation, as shown in FIG. 1, the desired output y^(k+ 
1) is given as: 

y(des)(k+ l)=y re /k+ l)+K Pe {y re /k)-y(k)}-K ie e i (k)-y ad 

{k) (6) 

The control input is given as: 


«(*) = {yref(k + 1) + K Pe {y ref (k ) - y(k)} - K^e^k) - y ad (k ) - ( 7 ) 

f A \y(k), y(k-l), ... , y(k -p y \ u(k- 1), ... ,u(k-p u )}}/ 

g{y(k),y(k- 1), ... , y{k-p y \ u(k- 1), ... ,u(k-p u )} 
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5 

To analyze the effect of this control input, we look at the 
modeling error, which is defined as the difference €(k+l) 
between the actual plant output and that predicted by the 
model: 

5 

e(k + 1) = y(k + 1) - (8) 

f A {y(k),y{k- 1), ... ,y(k- p y ;u(k- 1), ... , u(k - p u )} -g A 

{y{k),y(k- 1), ... , y{k - p y \ u{k - 1), ... , u(k - p u )} u(k) io 


Substituting the expression for the control input, given by Eq. 

(7), in Eq. (8) gives: 

15 

Ki e ei(k) (9) 

In terms of the definition of the tracking error, Eq. (9) can 
written as: 

20 

e(k+ l)+K Pe {y(k)-y re /k)}+K l e e l (k)=e(k+ 1 )-y ad {k) (10) 

Equation (10) represents a key equation of this approach. The 
left hand side of Eq. (10) is the desired second order error 
dynamics. The right hand side of Eq. (10) is the difference 
between the modeling error and adaptive augmentation signal 25 
input. Equation (10) indicates that, if the adaptive augmenta- 
tion signal can learn the modeling error and cancel this error, 
the error dynamics of this control loop will be restored to its 
desired nature. In other words, we will recapture the perfor- 
mance desired from this control loop. We, therefore, define 30 
the left hand side of Eq. (10) as the performance error, E(k), 
which is more realistically expressed as a vector of perfor- 
mance error components. 

Eik+l^eik+V+KpXfy+Kiee^k) (11) 35 

We can now form a Lyapunov function of the performance 
error as: 

L(k)^\E(k)\ 2 , (12) 4Q 

An update law now can now be devised for the adaptive 
augmentation input, y ad , that imposes monotonically 
decreasing behavior on this Lyapunov function. 
Parameterization and Update Laws for the Adaptive Augmen- 
tation. 45 

In this section, we investigate two questions. The first 
relates to the parameterization of the modeling error, and the 
second relates to the choices for designing stable update laws. 

1 . Parameterization for a Linear System: 

Consider a linear system of the form: 50 

x(k+l)=Ax(k)+Bu(k), (13) 

where x and u are vector components of the plant variables 
and the control inputs and A and B are system matrices. In a 
manner similar to that illustrated by Eqs. (3-5), the control 55 
input is computed as: 

u(k)=B '~ 1 {x re /k+ 1 )+K Pe e(k)+K le e : (k)-x ad (k)-A ~x{k) } (14) 

where A" and B~ are estimates of the system A and B matrices. 6Q 
If the system matrices (A, B) are known, adaptive augmenta- 
tion is not needed, and the control input is computed as: 

u(k)=B ~ 1 {x re /k+l)+K Pe e(k)+K le e : (k)-Ax(k)} (1 5) 

If these control inputs are to provide the same desired error 65 
dynamics, they must be equated, which gives the form of the 
idealized value of the augmentation signal x a ^k). 


B A 1 {x ref (k + 1) + Kp e e{k ) + K). e ei (k) - x ad (k) - Ax(k)}} = ( 16 ) 

B~ l {x ref {k + 1) + K Pe e(k) + K le e { (k) - Ax{k)} 


x A ad (k) = {(/ - B A B A - l )}{x ref (k + 1) 4- K Pe e(k) + K ie ey(k)} + (17) 

(B A B A ~ X A- A A )x(k) = 

{(/ -B A B A ~ 1 )(B A B A ~ l A — A A ){x ref (k + 1) + K Pe e(k) + K le e { (k)} 


2. Parameterization for a Non-Linear System Affine in Con- 
trol: 

Consider a non-linear system that is affine in control, and 
whose dynamics can be written as linear in parameters. 

x(k + 1 )=Wp/k)+Bu{k) ( 1 8) 

where Wyis the linear dynamic weight matrix, and the vector 
py corresponds to the linear and/or nonlinear functions of the 
system state. The control input is computed in a similar man- 
ner as: 

uWr'B-'^k+iy+K^ (19) 

where Wy " and B" are the corresponding estimates of the 
system matrices. By carrying out the analysis similar to the 
linear system case, the ideal augmentation signal can be com- 
puted to be: 


x ad (k) = {(I-B A B A ~ 1 )} 


(20) 


{(B A B A ~ l W f - W A f ){x ref (k + 1) + K Pe e(k) + K le e { (k)} 


Equations (17) and (20) imply that the ideal augmentation 
signal can be written as: 

p/*) (21) 

with the ideal weights, W*^, and the basis functions, p, as 
given in Eqs. (16) and (1 9). These are the same basis functions 
used in Rysdyk and Calise, ibid. Thus, we can parameterize a 
neural network in this form, and compute the ideal weights 
iteratively using an appropriate update algorithm. 

3. Update Laws for the Adaptive Augmentation. 

Having looked at the question of parameterization, we now 
construct a stable update law for the parameters W ad . Param- 
eterizing the adaptive augmentation signal in the form given 
by Eq. (20), and using the definition of the performance error 
as given in Eq. (1 1), Eq. (10) can be rewritten in vector form 
as: 


E(k+l)=e(k)-x a M, (22) 

Compared to Eq. (10), this is written for an error vector, E, 
corresponding to the general case of multiple control loops. 
Written in this form, the equation indicates that one estimates 
the vector modeling error, €(k) (for all loops) using the adap- 
tive augmentation signal x^k). The vector E(k+1) is a cor- 
responding error in the estimate. This error dynamics for the 
performance error E(k) corresponds to a system identification 
like problem. This opens up a host of approaches for doing 
this online system identification. In this work, we consider a 
normalized gradient update approach. 

4. Normalized Gradient Update. 

Let E z (k) correspond to the ith element of the vector per- 
formance error E(k). Let W*^ zz - represent the ith column 
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vector of the weight matrix W* flrf , which corresponds to the 
ideal weights that minimize the performance error vector 
components E z (k) to A*={6* l5 . . . , 6* z }. 

Similarly, let W ad>i represent the ith column vector of the 
current estimate of the ideal weight matrix. The update law 5 
for each of these column vectors of the weight matrix is given 
as: 


1)} (23) 

The parameter y (Eq. (12)) corresponds to the learning rate 
that lies in a range 

0<y=2 (24) 

Reference [11] proves that with this weight update law, the 
performance error, E(k). is monotonically decreasing for all i, 1 5 
Further, it is known that if the system experiences sufficient 
persistent excitation, the weights W arfl - approach the ideal 
weights W*^. 

5 . What Happens to Tracking Error? 

The final part of this analysis investigates the behavior of 20 
the system error e(k). This work provides an update only 
when modeling error is present, as opposed to presence of 
tracking error. However, tracking error is what is ultimately 
important. It is, therefore, appropriate to analyze the asymp- 
totic behavior of the tracking error given the behavior of the 25 
performance error. For simplicity, in this analysis we consider 
the case where the desired error dynamics is first order given 
as: 


(25) 30 

Let E.(k)<8. after time k, where 5 is some small positive 
scalar. This implies 

le 1 <£+l)-tf Pe e I .(A:)l< 6 , (26) 

From The Cauchy-Schwarz inequality, 35 

le^k+iyKp.e^l^le^k+l) I - IKpJeJk) I . (27) 

Equations (26) and (27) imply: 

40 

fo(k + 2)| < \K Pe \\ei(k + 1)| + 6, (28-1) 

<\K Pe \ 2 \e { (k)\+S{l + \Kp e \} 


\e i (k+n)\<\K Pe \ n \e l (k)\+b{l+\K Pe \+ . . . +\K Pe \ n ~ l} (28-n) 

Because IKp e l<l for stable error dynamics, ask— >oo, le.(k)l is 
bounded above as: 

\ ei (k)\<bl{\-\K Pe \} (29) 5Q 

Thus, if the performance error is bounded, Eq. (29) estab- 
lishes bounds on the tracking errors. A similar analysis can be 
carried out for second order error dynamics. The result sum- 
marizes that as long as the desired error dynamics (first or 
second order) is stable, the tracking error will be bounded 55 
above, given that the performance error is bounded. 

6. Application to Aircraft Control. 

The modeling error-driven performance-seeking adaptive 
control design was implemented for aircraft roll, pitch, and 
yaw rate control. The NASA Intelligent Flight Controller 60 
(IFC) design has been tested, and is currently undergoing 
various modifications for being flight-tested on the research 
F-15 aircraft. The IFC design has been implementing the 
adaptive control design as outlined by Rysdyk and Calise, 
ibid. For implementing the performance seeking adaptive 65 
augmentation, the requirement was that it needed to fit within 
the existing architecture. The main issue in the implementa- 


tion is that the baseline controller in the IFC architecture uses 
continuous-time aircraft dynamic inversion, whereas the pro- 
posed design has been outlined in the discrete-time. The 
equations outlined in the preceding sections have been for- 
mulated for a discrete-time model inversion. We realized, 
however, that after reducing the problem to the core error 
dynamics, the problems became identical. The error equation 
for the continuous -time implementation for a desired second- 
order error dynamics for a scalar error e is given as: 

(de/dt)+K Pe e+Kje(t')dt'=£.- U ad (k). (30) 

The error is defined in the same manner as the discrete case 
(e.g., q-q^e/). The modeling error, €, however corresponds to 
the difference in the acceleration, as predicted by the model, 
and the acceleration actually observed. Similarly, U ad repre- 
sents the augmentation acceleration command given by the 
adaptive block. If the left hand side of Eq. (30) is discretized 
while maintaining the continuous -time constants, the result- 
ing scalar discrete-time equation is given as: 


{e(k) + (K PE At-l)e(k - 1) + K le Atei(k - 1)}/A t = {s- U ad ). (31) 


Defining the left hand side of Eq. (31) as the modified per- 
formance error, E"(k), one obtains 

E\k)=e-U ad (k). (32) 

This modified performance error equation is identical to the 
discrete-time version given by Eq. (21). The adaptive aug- 
mentation acceleration signal U^k) can be parameterized in 
a similar manner, and the same update laws remain valid for 
the parameters of this augmentation signal for reducing E"(k). 
A zero value of this modified performance error restores the 
second order error dynamics (LHS of Eq. (31)) to zero, and 
thereby regains the desired performance from the control 
loops. Formulated in this manner, this adaptive approach fits 
within the existing IFC framework, and is considered as an 
alternate approach for flight testing. In the following discus- 
sion, we present some results of this implementation on the 
high fidelity model of the modified F-15 aircraft used at the 
NASA Dryden Flight Research Center. The adaptive control 
architecture is kept the same as in the original IFC design. 
This design has three loops, one for each of the pitch, yaw, and 
roll, respectively. Adaptive augmentation is provided to each 
loop. Kaneshige, and Burken., “Enhancements to a Neural 
Adaptive Flight Control System for a Modified F-15 Air- 
craft,” AIAA-2008-6986, give details on the implementation 
approach such as choice of the basis functions etc. The only 
difference is that the update law is given by Eq. (23). In this 
study, we look at two cases. In the first case, the right stabi- 
lator is locked at 4 degrees at t=10 sec into the flight experi- 
ment. In the second case, the canard multiplier is set at -1, 
again at t=10 sec into the flight experiment. The behavior of 
the aircraft and update algorithm is examined for the longi- 
tudinal and lateral pilot stick inputs given by FIGS. 2A-2B 
FIGS. 3A-3B through 7A-7D display the behavior of the 
aircraft and the neural net signal for right stabilator failure at 
t=10 sec. In FIGS. 3A-3C, the learning in the pitch channel is 
observable, when the aircraft actual pitch rate follows the 
desired pitch rate. FIGS. 4A-4C and 5A-5C display learning 
behavior in the roll and yaw channels, respectively. FIGS. 
6A-6C displays the behavior of the performance errors in 
each of the three axes as learning proceeds. These perfor- 
mance errors drive the updates in each of the three axe with 
performance errors becoming smaller as time increases. 
FIGS. 7A-7D display aircraft surface commands for this 
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maneuver. In a similar manner, FIGS. 8 A- 8 C, 9A-9C. 10A- 
10C and 11A-11C display the results for insertion of an 
incorrect canard multiplier. FIGS. 12A-12D display aircraft 
surface commands for the incorrect canard multiplier situa- 
tion. Note the good learning achieved in the pitch and roll 5 
channels. 

What is claimed is: 

1. A method for modeling error-driven adaptive control of 
an aircraft, the method comprising: 

providing a selected aircraft variable, y(k+l), at a time 10 
index having a value k+ 1 , as a matrix sum of Wy^k) and 
B u(k), where (3y(k) includes the at least one aircraft 
variable y(k) in a linear or nonlinear format, Wy is a 
matrix of selected aircraft variable weighting coeffi- 
cients, u(k) is a control variable vector for the at least one 1 5 
aircraft variable, and B is a matrix of control variable 
weighting coefficients that is not yet known; 
modeling an aircraft plant operation using a first neural 
network modeling mechanism, where the first neural 
network mechanism incorporates an assumption that the 20 
aircraft plant is operating within a normal range, without 
perturbations and without a tracking error vector e(k) 
that would cause the aircraft plant to experience an 
excursion outside a normal range of operation; 
providing a finite bound for the tracking error vector e(k) 25 
for operation of the aircraft within the normal range; 
when ( 1 ) at least one component of the tracking error 
vector e(k) is experiencing an excursion, determining if 
( 2 ) return of the at least one component of the tracking 
error vector e(k) toward a selected reference vector does 30 
not lie on or adjacent to a selected controller error char- 
acteristic; 

when the conditions (1 ) and (2) are satisfied for at least one 
value of the time index k, introducing at least one change 
in at least one parameter of the neural network modeling 35 
mechanism and modeling the aircraft plant operation 
according to a modified neural network mechanism with 
the at least one changed modeling parameter; and 
when the conditions (1) is satisfied and condition (2) is not 
satisfied, continuing to model the aircraft plant opera- 40 
tion using the first neural network mechanism, with little 
or no change in any modeling parameter of the first 
neural network mechanism. 

2. The method of claim 1, further comprising: 

expressing said tracking error e(k) as a difference 45 

e(k)=y^/k)-y(k) 

between a reference vector y re y(k) and said vector represent- 
ing said at least one aircraft variable y(k) at said time index k; 
expressing a change in said tracking error vector e(k) 50 
between said time index value k, and said time index 
value k +1 in a tracking error vector equation 

e (k+l)+K /?e e(k)+K 1<?ei (k)=0, 

where e : (k) represents an integrated tracking error up to and 55 
including said time index k, and K pe and K le are gain values; 
expressing a changed value in said aircraft variable y(k) as 
a sum 
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interpreting said control relationship as a linear or affine 
relationship 

y(k+l)=A y(k)+B u(k), 

where A is a matrix, with A as yet unknown, and said matrix 
B is invertible, with an inversion matrix B _1 ; 

providing estimates, A" and B", of the matrix A and said 
matrix B; 

expressing said control variable vector u(k) as a first sum 
u (k)=B _1 {y rey (k+ 1 )+K pe e(k)+K le e/k)}; 

expressing said control variable u(k) as a second sum 

u(k)=B^ 1 {y re /k + l) + K^e(k) + K le ei (k)-y^(k)- 
Ay(k)}, 

where y^k) is an adaptive augmentation vector that is 
included to encourage said aircraft to satisfy the tracking 
error vector equation: and 

expressing an alternative aircraft variable y ^(k) as a third 
sum 

/ B" B _1 )}{y re y(k+l)+K ??e e(k)+K le ei (k) 

}+(B~B -1 A — A~)y(k) ={(I — EfB -1 )(ITB -1 A — 

A") {y re Jk + 1 )+K pe e(k)+K le eftk)}. 

3. The method of claim 2, further comprising: 
expressing at least part of plant dynamics of said aircraft as 
a control vector relation, 

y(k+l)=f{y(k), y(k-l), . . . , y(k-p y ; u(k-l), . . . , u(k- 
Pj}++g(y(k), y(k-l), . . . , y(k-p^; u(k-l), . . . , 
u(k-pj} u(k), 

where f is a function of one or more of the variables y(k') 
(k ,= k, k-1, . . . , k-pj and one or more of the variables u(k'), g 
is a function of one or more of the variables y(k') and one or 
more of the variables u(k'), and f and g characterize the 
aircraft plant and are as yet unknown; 

inverting the control relation to express the control vector 
u(k) in a form needed to achieve a desired aircraft control 
dynamics as 

u(k)={ 7r< Jk+l) +K Pe { Yre/k) — y(k)} — K le e x (k) 

— f{y(k), y(k-l), . . . , y(k-p^; u(k-l), . . . , u(k- 
pj} Fg(y(k), y(k-l), . . . , y(k-p y ; u(k-l), . . . , 
u(k-pj}; 

providing an approximate model of the aircraft plant 
dynamics, with corresponding estimates, F and g", for 
the respective functions f and g; 
expressing an estimate u"(k) of a control vector that will 
achieve a desired aircraft control dynamics as 

uTk)={Y re /k+l)+K Pe { Yre /k) }— K x e e : (k)— y^(k) 

-f{y(k), y(k-l), . . . , y(k-p y ; u(k-l), . . . , u(k- 
Pj}FgTy(k), y(k-l), . . . , yik-p^; u(k-l), . . . , 
u(k- P J}; 

and 

providing an estimate e~(k) of said tracking error as said 
tracking error 

e (k+l)=y(k+l) — F (y(k), y(k-l), . . . , y(k-p y ; u(k-l), 

■ ■ ■ , u(k-p J}— g~{y(k), y(k-l), . . . , y(k-p y ; u(k- 
l),...,u(k- P JK(k). 


y(k+l)=y re /k+l) +K i7e {y r ^ )} .^ le e 1 (k); 





